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A significant, practical challenge, which 
arises in developing computationally effi- 
cient physical models tor use in computer 
simulations of microelectronic and op- 
toelectronic devices (for example, transis- 
tors in digital cellular phones and lasers 
in optical networks, respectively), is to rep- 
resent vast amounts of numerical data for 
transport properties in two or more dimen- 
sions in terms of closed form analytic 
expressions. In this paper, we present a 
general methodology to achieve the above 
goal for a class of numerical data in a 
bounded two-dimensional space. We then 
apply this methodology to obtain a closed- 
form analytic expression tor the minority 
electron mobilities at 300 K in p-type 



Gai-^Al^As as functions of the acceptor 
density Np, between lO" cm"' and lO"' 
cm"' and the mole fraction of AlAs x 
between 0.0 and 0.3. This methodology and 
its associated principles, strategies, re- 
gression analyses, and graphics are ex- 
pected to be applicable to other problems 
beyond the specific case of minority mobil- 
ities addressed in this paper. 
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1. Introduction and Motivation 



Researchers often first obtain theoretical and/or ex- 
perimental results in the form of multi-dimensional, 
numerical tables consisting of discrete data-points and 
then use families of traces to represent graphically their 
discrete data-points in a more easily understood manner. 
Such simplified graphical representations are a common 
recourse when several complex and competing physical 
mechanisms occur and when multidimensional, closed- 
form analytic expressions are not available. Incorporat- 
ing such discrete data tableaux into physical models for 
use in computer simulations is usually not satisfactory 
due to excessive computer time associated with interpo- 
lations between the discrete data points. Since such 
computational inefficiencies associated with "look-up" 
tables occur, industry is reluctant to incorporate them in 



semiconductor device simulators that run on engineer- 
ing workstations. This is particularly the case when nu- 
merical differences must be used to compute first and 
higher order derivatives. 

In this paper, we present a general strategy based on 
separable functions, melding functions, transformations, 
admissible non-linear methods, and regression analyses 
to obtain multi-dimensional, closed-form analytic ex- 
pressions from tables containing discrete data points. As 
an illustrative example, we apply this general strategy to 
a class of numerical data in a bounded two-dimensional 
space and show how to obtain a two-dimensional, 
closed-form analytic expression for the minority elec- 
tron mobilities in p-type Gai_^Al^As. The acceptor 
density A^a varies between 10"" cm"^ and 10™ cm"^ and 
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the mole fraction of AlAs x varies between 0.0 and 0.3. 
By so doing, we respond to the need for predictive 
computer simulations of devices that have fewer un- 
known or variational parameters [1]. 

Many device simulators for bipolar and field-effect 
transistors require a variety of physical models and asso- 
ciated input parameters to describe fully how carrier 
transport varies with carrier concentrations, ionized do- 
pant densities, alloy mole fractions, and temperature. 
This paper focuses on one such model, namely, the 
model for how the minority electron mobility varies 
with dopant density and mole fraction of AlAs at 300 K. 
Self-consistent numerical solutions to the quantum me- 
chanical, non-linear integral-differential equations for 
carrier transport in semiconductors result in discrete 
data-points that by themselves do not readily suggest 
closed-form analytic expressions for transport proper- 
ties, carrier densities of states, and band structure 
changes. Interpolating among the discrete data-points in 
"look-up" tables leads to discontinuities, particularly in 
numerical derivatives, and, as mentioned above, is com- 
putationally inefficient. 

The motivation for our performing the following anal- 
ysis is to derive a closed-form analytic expression that 
will result in more efficient computer simulations and 
improved physical insights on how the many scattering 
mechanisms, which influence carrier transport in 
ternary compound semiconductors and heterostructure 
devices, affect their electronic and optical behavior. Our 
data analysis, presented in the following sections, en- 
ables us to reduce the number of unknown physical 
parameters in numerical simulations that predict electri- 
cal and optical performance of devices such as bipolar 
transistors, solar cells, laser diodes, and light-emitting 
diodes. The latter are used to read digital versatile disks 
(DVD). 

The development of such improved closed-form mo- 
bility models directly impacts the design of microwave 
heteroj unction bipolar transistors used in the linear 
power amplifiers of digital cell phones. The design chal- 
lenges and goals are low noise and very linear, efficient 
power amplifiers. Such amplifiers enable longer talk 
times and improved adjacent-channel rejection in the 
dense channel packing that is necessary to maximize the 
capacity of communications systems. Designers of these 
amplifiers rely, in part, on device simulators to give 
them physical insights for selecting the optimum tech- 
nology and designing reliable, low cost devices. Using 
device simulators also provides sources of expert 
knowledge from others, shortens the time to market by 
reducing the number of experiments needed for design 
verification, and saves money. In order to increase the 
likelihood that computer simulations of heterostructure 
devices such as heteroj unction bipolar transistors 



(HBTs) will be predictive and thereby more useful for 
the future development of HBTs, it is essential to have 
accurate values for the minority mobilities of electrons 
in Gai-jAlvAs. This is especially true when one designs 
the bases of HBTs in mobile phones that use code divi- 
sion multiple access (CDMA) protocols. Such protocols 
require very linear power amplifiers for acceptable adja- 
cent-channel rejection. These amplifiers also should 
consume a minimum of energy. Device manufacturers 
may then have greater confidence when using more 
predictive simulators to design products. 

The markets for HBTs in mobile and lightwave com- 
munications systems are increasing [2-6]. For example, 
the worldwide number of new, cellular subscribers ex- 
ceeded the number of new subscribers for conventional 
wired telephones or fixed networks for the first time in 
1995 [7]. Market demands for telecommunications and 
optical storage systems that contain laser diodes (LDs) 
and light- emitting diodes (LEDs) offer many opportuni- 
ties for ternary III-V compound semiconductor HBTs. 
A common feature for most of these devices in such 
systems is that they contain active regions with high 
concentrations of dopants and carriers. Such high con- 
centrations alter the densities of states for the carriers 
and the band structure of the semiconductors. These in 
turn greatly affect how carriers move in semiconductors 
and modify substantially carrier transport parameters 
such as carrier mobilities, lifetimes of carriers, and ef- 
fective carrier concentrations. Accurate physical models 
for these transport parameters are needed in device sim- 
ulators so that they may be more predictive when man- 
ufacturers use them to design products. 



2. Data Table for Minority Electron 
Mobilities in p-type Gai-xAUAs 

This section summarizes the background details by 
which the data table for minority electron mobility val- 
ues was obtained. This data table serves as our starting 
point for deriving the closed-form analytic expression 
for the mobility model. 

The theoretical calculations in Ref. [8] include all the 
important scattering mechanisms for the low-field mo- 
bilities of electrons in heavily doped Gai_^Al^As: 
acoustic phonon, polar optic phonon, nonpolar optic 
phonon (holes only), piezoelectric, ionized impurity, 
carrier-carrier, plasmon scattering, and alloy scattering. 
The Boltzmann transport equation is solved by the vari- 
ational procedure outlined by Walukiewicz et al. [9]. 
This method avoids the use of the relaxation-time ap- 
proximation that is invalid for mechanisms that involve 
energy transfers comparable to or greater than k-sX, 
where k-n is the Boltzmann constant, and T is thermody- 
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namic temperature. Also, Matthiesson's rule is not used 
because it is not valid for Gai-^ Al;i As, as shown in Ref. 
[9] . The scattering rates are summed prior to the varia- 
tional solution. The results are highly accurate calcula- 
tions of the minority mobilities for p-type Gai-jAlvAs 
at 300 K as functions of the dopant density and the mole 
fraction of AlAs. The calculated results agree well with 
the rather limited experimental data for Gai-jAlvAs 
mobilities [10-12]. 

The calculations were done for a full factorial design 
consisting of 21 discrete values of acceptor density A^a 
between 10"" cm"^ and 10^" cm"^ and seven discrete 
values of mole fraction x between 0.0 and 0.30, namely, 
x = 0.0, 0.05, 0.10, 0.15, 0.20, 0.25, and 0.30 [denoted 
also by X = 0.00 (0.05) 0.30, a general notation that we 
use latter on in this paper], to yield a total of 147 data 
points. We use the notation that Xi = 0.00, X2 = 0.05, ... 
Xy = 0.30, respectively. The self-consistent, numerical 
solutions to the quantum mechanical, non-linear inte- 
gral-differential equations for carrier transport in semi- 
conductors are given in Table 1 as a two-dimensional 
array of discrete data-points. This data representation, as 
opposed to a functional representation, was necessary 
because the several competing scattering mechanisms, 
listed above and considered all together, do not readily 



yield any acceptable theoretical closed-form analytic ex- 
pression. Over 90 h of NIST Cray YMP supercomputer' 
time were used to obtain these 147 data points. 

The 147 data points presented in Table 1 are repre- 
sented graphically in Fig. 1 as a family of seven traces 
corresponding to the seven mole fraction values x = 0.00 
(0.05) 0.30, respectively. The fixed increment of 
x,_i — X, = 0.05 for all / and a subsequent fortuitous 
response surface in the mole-fraction variable will be 
advantageously employed later to simplify the fitting 
process. 

We thus have the task of finding a closed-form two- 
dimensional analytic function g for the minority elec- 
tron mobility fx^ in p-type Gai_^Al;iAs such that 
yu-c = g(NA, x). For simplicity we may alternatively find 
a function / for the dimensionless normalized mobility 
Y = jJbJiJbrcf such that 



y=/(^A,x) 



(1) 



where g = ju,rcf/and /^rct = 1000 cmV(V-s). To achieve 
an acceptable analytic fit of Y to A^a and x for use in 
semiconductor device simulators that run on engineer- 
ing workstations we require a relative residual standard 
deviation for Y that is less than 2 %. 



Table 1. Two-dimensional array of data points from tlieoretical calculations of tlie normalized minority electron mobility" 



Wa 








X 








10'^ cm-' 


0.00 


0.05 


0.10 


0.15 


0.20 


0.25 


0.30 


1.00 


5.742 


5.136 


4.552 


4.003 


3.500 


3.051 


2.658 


2.00 


5.106 


4.584 


4.081 


3.606 


3.171 


2.779 


2.434 


3.00 


4.718 


4.244 


3.788 


3.357 


2.961 


2.604 


2.288 


5.00 


4.225 


3.810 


3.411 


3.034 


2.686 


2.371 


2.092 


7.00 


3.903 


3.524 


3.161 


2.818 


2.501 


2.214 


1.959 


I.OOXIO' 


3.565 


3.224 


2.897 


2.588 


2.303 


2.045 


1.814 


2.00X10' 


2.926 


2.654 


2.394 


2.148 


1.921 


1.715 


1.530 


3.00X10' 


2.569 


2.334 


2.110 


1.899 


1.703 


1.525 


1.366 


5.00X10' 


2.145 


1.953 


1.771 


1.599 


1.441 


1.296 


1.165 


7.00X10' 


1.889 


1.721 


1.564 


1.416 


1.279 


1.154 


1.041 


1.00X10' 


1.643 


1.499 


1.364 


1.239 


1.123 


1.016 


0.9198 


2.00X10' 


1.258 


1.150 


1.050 


0.9575 


0.8725 


0.7948 


0.7241 


3.00X10' 


1.101 


1 .006 


0.9191 


0.8397 


0.7670 


0.7007 


0.6403 


5.00X10' 


0.9936 


0.9046 


0.8249 


0.7530 


0.6880 


0.6292 


0.5760 


7.00X10' 


0.9901 


0.8971 


0.8152 


0.7420 


0.6764 


0.6175 


0.5647 


I.OOXIO' 


1.068 


0.9601 


0.8660 


0.7832 


0.7097 


0.6449 


0.5873 


2.00X10' 


1.539 


1.372 


1.227 


1.098 


0.9863 


0.8881 


0.8013 


3.00X10' 


1.964 


1.746 


1.556 


1.390 


1.243 


1.115 


1.002 


5.00X10' 


2.489 


2.212 


1.968 


1.753 


1.563 


1.397 


1.250 


7.00X10' 


2.692 


2.392 


2.125 


1.892 


1.687 


1.506 


1.347 


1.00X10" 


2.720 


2.418 


2.152 


1.918 


1.711 


1.528 


1.368 



' Y = /ic(p-type; Wa, A:)//j,rei for p-type Gai-,Al^As, where the acceptor density is A'a, the mole fraction of AlAs is x, and /iref = 1000 cm /(V-s) [8]. 

' Certain commercial equipment instruments or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available tor the purpose. 
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Fig. 1. Minority electron mobility in p-type Gai-<Al<As for several values of mole fraction in 
the range 0.0 < j: < 0.3. These curves include the effects of the electrons scattering off hole 
plasmons and of deducting the density of holes with energies below the Fermi energy from the 
electron-hole scattering process [8]. The plasmon cut-off factor (PCF) q^r^ is 1 where q^ and r, 
are the cut-off wave number and the screening radius, respectively. The mobilities have been 
normalized to ij.„f = 1000 cmV(V-s). The dimensionless acceptor density is 7V„ = N^ cm^. 



The development of such a function /would represent 
a significant increase in computational efficiency by 
about a factor of 5 and gives mobility models for use in 
commercial semiconductor device simulators that are in 
closer agreement with measurements. The combination 
of the existing NIST supercomputer-generated mobility 
data and the derived two-dimensional analytic function 
/ will lead to computer simulators that are at once both 
more parsimonious (have fewer unknown or tuning- vari- 
ational parameters) and more accurate (offer improved 
predictability). 



3. Data Analysis for Minority Electron 
Mobilities in p-type Gai-^ALAs 

In the following sections we show that using a combi- 
nation of separable functions, melding functions, trans- 
formations on the discrete data points in Fig. 1, and 
non-linear regression analyses leads to a single two-di- 
mensional, closed-form analytic expression for the mi- 
nority electron mobilities at 300 K in p-type 
Gai_j;Alj^As as functions of the mole fraction of AlAs x 
between 0.0 and 0.3 and the acceptor density A^a be- 
tween 10"" cm"' and 10™ cm"^ Throughout our analyses 
we rely substantially on graphics and keep the number of 



fitting coefficients to a minimum, subject to the con- 
straint that the residual standard deviation Sr^, (Y) as 
defined by Eq. (2) satisfies 5'rcs(I') ^ 0.02. The residual 
standard deviation is a measure of the "average" error in 
a fitted model and thereby is a metric for assessing the 
quality of the fit, with a smaller Src^iY) indicating a 
better fit. The residual standard deviation for a model 
Y=f(Nj,,x) is 



SUY) = \ [^{Y, - Y.f/in - p)] 



(2) 



where Yj are the observed data values, the Yj are the 
predicted values from the fitted model, n is the total 
number of data points (here n = 147) and p is the total 
number of parameters to be fitted in the model (here 
p = 18). The residual standard deviation for a model, Eq. 
(2), with p parameters differs from the standard devia- 
tion for a set of data points So(Y). The standard deviation 
for a set of data points is given by 



Sn(Y)=J[ZiYj-Yf/{n- 1)] 



j=i 



where Y=^Yj/n is the arithmetic average. 
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We use the NIST-developed DATAPLOT [13] soft- 
ware for both the exploratory graphics and for the exten- 
sive non-linear statistical analyses. Also, for those cases 
in which the residual standard deviations from analyses 
based on different functional forms are quantitatively 
similar, we select the functional form that will minimize 
the computer time when the closed-form analytical 
function is used in commercial simulators and select 
procedures that have a minimum of fitting parameters. 

Our general strategy is based on separable and meld- 
ing functions and on transformations of the response 
function Y that give near-linear separable functions as 
described below. We want to obtain the function 
Y=f(X,x) in the two-dimensional continuum space 
spanned by X = log[A'A/(10"' cm"^)] and x. This 
bounded two-dimensional continuum is given in Fig. 1 
with < X < 4 and 0.00 < x < 0.30. As with fitted 
functions, extreme caution must be exercised in extrapo- 
lating beyond these X and x limits. 

3.1 Separable Functions 

Consider the discrete two-dimensional space given by 
the 147 data points in Fig. 1. Let F, denote the ith data 
trace in Fig. 1 with / = 1(1)7 corresponding to mole 
fraction x = 0.00(0.05)0.30. //the 7, are related to F, via 
expressions of the form 

Y, = H,iYO 

and if Yi is related to X by 

Y,=BiX) 

then Y, = Hi[BiX)]. 

In practice, //, and B are least-squares, best-fit func- 
tions relating F, to Yi and Yi to X, respectively. We call 
the functions F and B separable functions. In addition, 
if the F were to be linearly related to Fi, that is, if 

Y, = H,iY,) = h + SiYu 

where /, and 5, denote the j'th intercept and the / th slope, 
respectively, then 

Yi = Hi[B{Xy[=h + S,B{X). 

In this case, F and B are called linearly separable func- 
tions. 

We now determine the extent to which the raw data in 
Fig. 1 may be represented by separable and by linearly 
separable functions. Without any additional information 
besides that contained in Fig. 1, we could select any one 
of the seven traces as a base trace and examine how the 



other six traces relate to the base trace. However, since 
we know that trace Fi for GaAs (the x = 0.0 trace in Fig. 
1) has been verified experimentally by three separate 
groups [10-12], we can use with no loss in generality, 
trace Fi as the reference or base trace. 

We now discuss the strategy for determining the na- 
ture of the relationship //, between the trace F and the 
base trace Fi. We tentatively hypothesize that the trace 
F, is linearly related to Fi. We would thus have the 12 
parameter representation for / = 2(1)7, 



F = H,(F,) = /, + 5,F,^ 



(3) 



with six values of intercepts /, and six values of slopes 
Si. The graphical analog for Eq. (3) is a fan-shaped 
collection of traces in the F vs Fi plot. 

The first simplification for Eq. (3) would be for all six 
slopes Si to be identically equal to a constant slope S; 
thus, yielding the seven parameter representation 



Yi = h + SY,. 



(4) 



Graphically, Eq. (4) results in a F versus Fi plot consist- 
ing of seven parallel lines. 

The second simplification would be for the six inter- 
cepts to be a constant multiple of / as in 
Ii = I+{i— l)D, where / and D are constants, thus 
yielding the three parameter representation 



F = /+(('- l)D + SYi. 



(5) 



Graphically, Eq. (5) results in a F vs Fi plot consisting 
of seven parallel lines having constant spacing. 

Examination of Fig. 2 shows that the seven traces of 
F vs Fi for the mobility data are not perfectly parallel 
and do not exhibit fixed spacing. This indicates that the 
simple three parameter representation is not adequate. 
We note from Fig. 2 that the family of seven traces do 
exhibit rough linearity with different intercepts /, and 
with markedly different slopes Si, thereby tenuously 
justifying our 12 parameter linear representation. At this 
point a series of six linear fits could be carried out from 
Fig. 2, thus yielding estimates for the 12 parameters (six 
intercepts and six slopes). 

We did not do this for the following four reasons: 

1. Non-linearity: From Fig. 2, we note that because the 
traces exhibit subtle but progressively increasing cur- 
vature for / = 2(1)7, the linearity assumption is vio- 
lated. 

2. Parsimony: Parsimony dictates that the number of 
parameters be as small as possible and it is preferable 
to reduce the 12 parameters to a smaller number. 

3. Stability: Estimating fewer parameters is usually 
computationally more stable. 
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Fig 2. Plots of YiiX) vs ¥,{X) for ; = 1(1)7. 



4. It may be possible to find a representation for the 
p-type minority electron mobilities that is more fun- 
damental and hence may be applicable to other semi- 
conductor modeling cases. 
For the reasons cited above, we set aside the 12 

parameter representation and focus on a modified (and 

ideal) three parameter representation. Clearly, from Fig. 

2, this cannot be achieved for the original normalized 

values of Y, because 



Yi^ I+(i- 1)^-1-571. 



(6) 



Our strategy is then to seek in Sec. 3.2 the transforma- 
tion U of the mobility data such that 



Z, = In Y, = lim[(y/ 



1)/A]. 



(10) 



Figures analogous to Fig. 1 for each of the five values of 
A show that when A = 0, the seven traces of Z vs X are 
the closest to being nearly parallel and equally spaced as 
shown in Fig. 3. Figure 4 then gives the seven traces of 
Z, vs Zi. The next step is to plot in Figs. 5 and 6 the 
intercepts /, and the slopes 5, for the transformed traces 
Z,(Xq) = In Yi(Xo) at X = Xo = vs /. Polynomial fits for 
the intercepts and slopes give expressions of the form 



/(X;) = flo + fll Xj + ^2 Xi^ + flj X;^ + ... + a„ Xj' 



(11) 



and 



UiY,) = I+(i- l)D + SU(Yi). 



(7) 



We assign the notation Z, to the transformed mobility 
data U{Yi) and obtain 



Z, = U(Yi) = I+(i- l)D + SZi. 



3.2 Transformation of Mobility 



(8) 



To determine the appropriate transformation U we 
consider the following extended power-transformation 
family: 



Z = (F/ - 1)/A 



(9) 



for A = — 1.0(0.5)1.0. When A = the transformation is 
the natural logarithm, namely 



S(Xi) = So + Si Xi + S2 Xi^ + S3 Xi^ + . 



.. + S„ Xi 



(12) 



where jc, = 0.05 (/ — 1) and / = 1(1)7. The condition of 
near-constant spacing for /, means that lao + aiXil 
should be much greater than 1 02 xi^ + 0^X1^ + ...+ a„ x,"l; 
and the condition of near-constancy for 5, means that I Sa\ 
should be much greater than \si x, + S2Xi^ + s^Xi^ + ... + 
s„ x,"l. The residual standard deviations 5'rcs(/) and SrcsiS) 
of polynomial fits to I(x) and S(x) are calculated for 
increasing orders n of the polynomials until acceptably 
small values of Sr^s occur. For the data in Fig. 1, cubic 
fits to both I{x) and S(x) give 5„,(/) = 0.000 11 and 
SmsiS) = 0.000 04 in the logarithmically transformed 
space. These values of 5'rcs(/) = 0.000 11 and 5'rcs(5') = 
0.000 04 are sufficiently small to achieve the overall 
goal of SUY) ^ 0.02. 
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Fig. 3. Plots of the natural logarithm of the normalized minority electron mobility in p-type 
Gai-, Al, As for several values of mole fraction in the range 0.0 < j: < 0.3. These curves include 
the effects of the electrons scattering off hole plasmons and of deducting the density of holes 
with energies below the Fermi energy from the electron-hole scattering process [8]. The plasmon 
cut-off factor (PCF) q^r^ is 1 where q^ and r, are the cut-off wave number and the screening 
radius respectively. The mobilities have been normalized to /J,,ei = 1000 cmV(V-s). The dimen- 
sionless acceptor density is N„ = A'a cm^ 
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Fig. 4. Plots of the transformed function Zi(X) vs Z,(X) for / = 1(1)7. 
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Fig. 5. Plot of the intercepts /, for the transformed curves Z|(Xo) at X = Xq vs j = 1(1)7. The mole 
fraction x is given by ;c = 0.05 (i — 1). 
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Fig. 6. Plot of the slopes S, for the transformed curves Z, (Xo) at X = Xo = vs ;' = 1(1)7. The mole 
fraction x is given by a: = 0.05(/ — 1). 
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In summary, the general strategy is to replace the 
diff icuh task of finding one function / with three easier 
tasks of finding a base function B (X), an intercept func- 
tion lix), and a slope function S{x) that is weakly depen- 
dent on X and for which / is the surrogate for the mole 
fraction variable x = 0.05(i — 1). 

3.3 Melding Functions 

The task of finding the reference or base function 
B(X), which is essentially the minority electron mobil- 
ity for GaAs, involves the following series representation 
for B(Xy. 



L(X, Xo,(To) = l/{ 1 + exp[-(X - Xo)/a-o] j 



(17) 



B(X) = ^w,(X)g,(X), 



(13) 



where the weighting functions Wi{X) satisfy the sum 

N 

rule 2 Wi(X) = 1. The function B{X) is represented by 

a melding function. A melding function in this case is a 
global A? -region additive function that joins the A^ local 
functions gi{X) to give a smooth function by the use of 
the weighting functions Wi{X) that are bounded between 
1 and over the global region of X. For the illustrative 
example here X is bounded between and 4. 

Figure 1 suggests that it is necessary to consider only 
two regions namely N =2. That is, we write B{X) as 



BiX) = wiX) g,(X) + (1 - wiX)) g2iX). 



(14) 



Commonly used functions for the local functions gi(X) 
are polynomials of nth order, rational functions, spline 
functions, other flexible functions, and melding func- 
tions themselves. We use cubic polynomials for the 
functions gi{X), namely. 



and 



gi(X) = Cio + CnX+CuX' + CuX' (15) 



gliX) = C20 + CnX+ C22 X' + C23 X\ (16) 



because they, as for the above intercepts and slopes, give 
acceptably small Sr^^igi) and Sr^^igz) values. 

There are initially many choices for the weighting 
functions w(X). The requirement that 5'rcs(}') ^ 0.02 
eliminates many of them. The two weighting functions 
that give the smallest Si^^{Y) values for the data in Fig. 
1 are the normal cumulative distribution function (CDF) 
and the logistic CDF (LCDF). The S^Y) values for the 
normal CDF and the LCDF are almost the same. Since 
the LCDF requires less computer time to evaluate, we 
select it for the remaining analyses. The LCDF is given 
by 



where Xo and (Tq are location and scale parameters of the 
logistic distribution. 



4. Final Results — Closed-Form Analytic 
Function 

The function B (X) is a global function in InY space 
composed of the two local functions gi{X) and giiX) 
that have the same value at the region boundary X = Xb. 
Figure 1 suggests that Xb should be near Xb = 2.75. The 
next step is to obtain cubic non-linear fits for gi(X) over 
the region ^ X < Xb and for giiX) over the region 
Xb < X < 4.0. 

The regression analyses involves three major steps: 

1. Determine the non-linear base function B{X) that 
relates the logarithm of the reference mobility trace 
Zi directly to X. 

2. Determine the intercept function /, and the slope 
function 5, that relate the other logarithms of the 
mobility traces Z, to Zi. In effect this allows the 
mole-fraction variable x to enter the fit by means of 
the relation x = 0.5(/ — 1). At this point, we have 
adequate non-linear fits relating the Z, to the doping 
density X and the mole-fraction x via the expressions 
Zi = H,[B(X)]=Ii + Si[B(X)]. 

3. Use the estimated fitting coefficients from Step 1 
above for the base function fitting [namely Cy, C2j, Xo 
and (Jo, wherej = 0(1)3] as initial values for the final 
non-linear fit relating the mobilities 7, to X and x by 
means of the function 



Y, = exp[Z] = exp{Hi[B(X)]} = exp[/, + SiB(X)]. (18) 

For this final fit, only the 10 base function coeffi- 
cients Cy, C2j, Xa and cto are varied. There is no need to 
vary again the intercept and slope coefficients Oj and Sj 
from Step 2. 

Combining all of the numerical and statistical proce- 
dures described in Steps 1 to 3 above, we obtain a 
two-dimensional closed-form analytic expression for 
electron mobilities in p-type Gai_^Al^ As as functions of 
acceptor density A^a and mole fraction x of AlAs, 
namely. 



/j,e(p-type; A'a, x) = fji,^fexp{I(x) 
+ S(x)[wL(X)g,(X) + (1 - wdX))g2(X)]] 



(19) 



where A^a = 10^ X 10"" cm""*, where the seven traces 
containing 147 data points in Fig. 1 are represented by 
the four cubic polynomial functions /(x), S{x), gi{X), 
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and giiX), and the exponential type function 
Wl(X) =1 - L(X, Xq, (To) for a total of 18 fitting coeffi- 
cients and 5'rcs(i') = 0.018. Table 2 contains the final 
values for the 18 fitting coefficients and other data that 
will be discussed in the next paragraph. 

First and foremost, the Srcs{Y) = 0.018 is sufficiently 
small for the purposes of simulating the electronic and 
optoelectronic behavior of semiconductor devices, and 
verifies the adequacy of the fitted function, Eq. (19), and 
its 18 coefficients given in Table 2. The value of 
Sics{Y) = 0.018 corresponds to a relative residual stan- 
dard deviation of about 10 % which is more than 
adequate for most microelectronic and optoelectronic 
computer simulations. Second, the condition of near- 
constancy of spacing of the intercepts /, is met because 
for the worst case of x = 0.30, I Aq + fli JcUo.3 = 0.593 is 
much greater than \a2X^ + a^ x\=().3 = 0.038. And third, 
the condition of near-constancy of the slopes 5, is met 
because for the worst case of x = 0.30 I i'ol = 1 .0 is much 
greater than \six + Szx'^ + s^x\^a3 = 0.\99. All 10 of 
the coefficients in B{X) are significant because the ab- 
solute values of their estimated standard deviation ratios 
greatly exceed the nominal value of 2. The estimated 
standard deviation ratio is the estimated value of a fitting 
parameter divided by its estimated standard deviation. It 



gives the number of standard deviations by which a 
fitting parameter differs from zero. The nominal value 
of 2 is based on the 97.5 % point of the normal (Gaus- 
sian) distribution. In practice, coefficients with esti- 
mated standard deviation ratios less than two are consid- 
ered to be statistically insignificant. Because the 
estimated standard deviation ratio of the intercept fitting 
parameter Uq is less than 2, Uq could be omitted from the 
final model with little effect. In fact, setting Ua = 0.0 
gives the same value Sr^siY) = 0.018 to within 10 ''. The 
collective significance for the 17 remaining fitting 
parameters implies that the fitted function is parsimo- 
nious in its form and cannot be "trimmed" to a simpler 
form with fewer parameters or coefficients. For a fur- 
ther discussion of the above, we refer the reader to 
Draper and Smit [14]. 

Conventional non-linear fitting procedures usually do 
not contain regression analyses, transformations, and 
separable functions. As a consequence, they typically 
require more than 40 fitting coefficients for the data in 
Fig. 1, yield Sr^siY) values for the data in Fig. 1 that are 
greater than 0.07. They therefore give analytic fits to the 
minority electron mobilities that may be of questionable 
value for use in device simulators. 



Table 2. The 18 final fitting parameters for the minority electron mobility from Eq. (19)° 



Intercept and slope 
fitting parameters 

a, 
ai 
ai 
So 

Sl 
Sl 
Sl 



Estimated 
value 

-0.000 05 
-1.976 36 

0.524 177 
-0.339 191 

0.999 986 
-0.061 328 
-1.530 36 

1.695 421 



Estimated standard 
deviation 

0.1052 X 10"' 
0.3323 X 10"^ 
0.2716 X 10"' 
0.5943 X 10"' 
0.3411 X 10"* 
0.1077 X 10"^ 
0.8808 X 10"^ 
0.1927 X 10"' 



Ratio 



-0.48 

-0.59 X lO-'^ 
19.0 

-5.7 
2.9 X 
-57.0 
-1.7 X 
88.0 



10*' 



10* 



Base-reference function 
fitting parameters 

Cll 

Cn 
Co 

C20 
Cll 
Cll 
Cli 

X„ 
era 



1.741 24 
-0.350 382 
-0.076 087 
-0.038 110 
68.654 8 
-50.516 9 
12.841 4 
-1.106 75 
3.242 62 
0.259 828 



0.001 598 
0.009 642 
0.014 74 
0.007 677 
0.103 7 
0.216 
0.002 999 
0.011 30 
0.042 98 
0.009 315 



1100. 

-36. 

-5.2 

-5.0 

660. 

-230. 

4300. 

-98. 

75. 

28. 



" /j,e(p-type; Nn, x) = iJ.„fexp{I(x) + S(x)[wi^(X)g,(X) + (1 — WL(X))gi(X)]], where I(x) = ua + a\ x + Ui x^ + Ui x^, S(x) = .«o + .S'l x + Si x^ + .«3 x', 
g\(X) = Cio + t| 1 X + Cll X^ + C|3 X^, gi(X) = Cia + CuX + Cn X^ + Cn X^, and WiiX) = 1/(1+ exp[(X — Xo)/a-o] } where the normalized-dimension- 
less acceptor density is X = logio (A'a/IO"' cm"'), the mole fraction of AlAs is x, and /iref= 1000 cmV(V-s). All of the fitting parameters are 
dimensionless. The ratio is the estimated value divided by its estimated standard deviation. 
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The above closed-form fit of the data in Fig. 1 is of 
such precision that the fitted curves lie almost within the 
Hne-widths of the theoretical traces. The analytic ex- 
pression in Eq. (19) now enables quantum mechanically 
based results, which required tens of hours of supercom- 
puter time, to be readily and efficiently incorporated 
into commercial workstation-based simulations of 
HBTs. 

The analytic fit in Eq. (19) is valid only within the 
ranges ^ X < 4 and ^ x ^ 0.30, and must not be 
used beyond this bounded two-dimensional space in 
which it is derived. Also, combining Eq. (19) with other 
transport models for mobilities, bandgaps, and effective 
intrinsic carrier concentrations that are derived from the 
interpretation of electrical measurements on the devices 
themselves may lead to incorrect descriptions of the 
electrical and optical behavior. 



5. Engineering Significance 

Using the above Eq. (19) and applying additional 
results from recent calculations [8,15] to microwave 
HBTs [16] for linear power amplifiers may suggest dif- 
ferent design strategies to optimize HBT performance. 
The calculated changes in carrier densities of states 
(DOS), band edges, band offsets, effective carrier con- 
centrations «ic, and carrier mobilities due to high dopant 
and carrier concentration effects in Gai-^jAl.As are 
given in Refs. [8] and [15] at 300 K for mole fractions 
X of AlAs between 0.0 and 0.3, for donor densities A^d 
between 10"" cm "* and 10''' cm ^ and for acceptor den- 
sities Na between 10"" cm"^ and 10^° cm ^ Only one 
quantum mechanical theory is used to treat both sides of 
the Mott transition in these calculations. They give, with 
no fitting parameters to experimental measurements, an 
internally self-consistent description of carrier transport 
in Gai_;iAl;iAs/GaAs heterostructures for lasers, light 
emitting diodes, digital devices, and microwave devices. 
The predicted values for the distorted DOS, band edges, 
band offsets, «!„ and majority and minority mobilities 
differ significantly from those values found in many 
simulations of Gai_vAlvAs/GaAs heterostructures. 
Many simulators set njrii = 1 in Gai-^i Al^jAs for all A^d 
or A'a, approximate /Xc(p-type; A^a) with /.tc(n-type; 
A^D = A^a), and assert that all mobilities are monotoni- 
cally decreasing functions of the dopant density. How- 
ever, Fig. 1 shows that a relative minimum exists for 
/J-e(p-type; A^a) and suggests that a different design strat- 
egy could be significant for linear HBT amplifiers in 
digital cellular phones. Because a relative minimum in 
the minority electron mobility as a function of the ac- 
ceptor density exists, we have identified additional de- 
sign considerations for HBT power amplifiers that 
would have not otherwise been known. The above rela- 



tive minimum arises from dependencies of several com- 
peting scattering mechanisms on both the dopant and 
carrier densities. A reduction in the scattering of minor- 
ity electrons off hole plasmons and the removal of ma- 
jority hole carriers from minority carrier-majority car- 
rier scattering due to the Pauli exclusion principle 
accounts for the relative minimum in the decade of 
10" cm I 

If other parameters remain essentially the same and 
A^A increases from 6 X 10"* cm"^ to 6 X 10''' cm"^^ then 
the following occurs: 

1 . the minority electron mobility increases by a factor 
of 2.5 [8], 

2. the base transit time decreases by about a factor of 
2.5, and 

3. the base resistivity decreases by about a factor of 10 
[17]. 

Combining these last three results into expressions from 
compact models for microwave HBTs predicts increases 
in operating frequencies of about 40 % and in figures of 
merits (maximum frequencies at unity gain) of about 
300 %. These estimates are considered to be upper lim- 
its because more rigorous simulations depend on both 
processing and operating parameters whose choices are 
determined by the application. 



6, Conclusions 

We have constructed a two-dimensional closed-form 
analytic function for the minority electron mobilities in 
Gai_jAljAs at 300 K that is a function of the mole 
fraction and acceptor doping density. All of the impor- 
tant scattering mechanisms including carrier-carrier and 
plasmon scattering are considered. The minority elec- 
tron mobility from first-principles quantum mechanical 
calculations shows an interesting structure at high densi- 
ties due to the reduction in plasmon scattering and the 
Pauli exclusion principle. The results are important for 
device modeling because of the need to have accurate 
values for minority mobilities which in turn allow im- 
proved design of HBTs for microelectronic and op- 
toelectronic applications, for example, digital cellular 
phones and modulators in optical communications sys- 
tems, respectively. 

The general modeling approach given here may be 
used to: 

1. generate close-form analytic expressions 

2. eliminate interpolator look-up tables and 

3. enable more efficient computer simulations 

for many microelectronics and optoelectronic applica- 
tions such as obtaining expressions for carrier transport 
properties and extracting two-dimensional doping pro- 
files from scanning capacitance microscopy measure- 
ments. 
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